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Abstract 

I ■ B^B^ flavor oscillations have been studied in 20.7 fb~ of e^e" annihilation data collected in 1999 

Qh. and 2000 with the EABAR detector at center-of-mass energies near the T{AS) resonance. The event 

r~| ! sample consists of one B^ meson fully reconstructed in a hadronic decay mode, while the flavor of 

the recoiling B^ in the event is determined with a tagging algorithm that exploits the correlation 

- between the flavor of the heavy quark and the charges of its decay products. By fitting the time 

r\ \ development of the observed mixed and unmixed final states, the B^B^ oscillation frequency, Am^, 

C^ ■ is determined to be 0.519 ± 0.020 ± 0.016 ?ips~^ 
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1 Introduction 

In the Standard Model, B^B^ |j^ mixing occurs through second-order weak diagrams involving 
the exchange of up-type quarks, with the top quark contributing the dominant amplitude. A 
measurement of Attz^, the difference between the mass eigenstates of the B^~B^ system, is therefore 
sensitive to the magnitude of the Cabibbo-Kobayashi-Maskawa matrix H element Vtd- At present 
the sensitivity to Vtd is not limited by experimental precision on Arrid, but by other uncertainties 
in the calculation, in particular the quantity f^Bs, where Jb is the B^ decay constant, and Bb is 
the so-called bag factor, representing the strong interaction matrix elements. 

The phenomenon of particle-anti-particle mixing in the neutral B meson system was first ob- 
served almost fifteen years ago [^. The oscillation frequency has been extensively studied with 
both time-integrated and time-dependent techniques |^. 

In this paper we present a measurement of time-dependent mixing based on a sample of 20.7 fb^ 
of data recorded at the T(4S) resonance with the BABAR detector at the Stanford Linear Accelerator 
Center. At the PEP-II asymmetric-energy e~^e~ collider, resonant production of the T{4S) provides 
a copious source of B^B^ pairs moving along the beam axis {z direction) with a Lorentz boost of 
P'j = 0.56. The typical separation between the two neutral B decay vertices is (|A2;|) « (3^ctbo = 
260 /xm, where r^o = 1.548 ± 0.032 ps is the 5° lifetime [§. 

2 Analysis method 

The B^B^ mixing probability is a function of Am^ and the time difference between the B decays. 
At ~ Az/p-fc: 

Prob{B°B^ -^ B^B°, WB^) cx e~l^*l/^so (i _ cos Am^At) 
Prob{B°B° -^ B^B^) cx q-\^^\/^bo (i + cos ArrirfAt) 

resulting in a time-dependent probability to observe mixed, B^B^ and B^B^ , or unmixed, B^B^, 
pairs produced in T{4:S) decay. The effect can be measured by reconstructing one i? in a flavor 
eigenstate, referred to as B^ec, while using the remaining particles from the decay of the other B, 
referred to as Btag, to identify, or "tag", its flavor. The charges of identified leptons and kaons are 
the primary indicators, but other information in the event can also be used to identify the flavor 
of Btag, resulting in a total of four non-overlapping tagging categories. The tagging algorithm 
used in this analysis is identical to that employed for CP violation studies, in which one B is fully 
reconstructed in a CP eigenstate [^. 

If the flavor tagging were perfect, the asymmetry as a function of At 



-'^mixing l,^^ j 



^^ unmixed ^^mixed 
-^^ unmixed ~r -'■^mixed 



would be a cosine with unit amplitude. However, the tagging algorithm incorrectly identifies the 
flavor of Btag with a probability Wi for the i*'^ tagging category. This mistag rate reduces the 
amplitude of the oscillation by a factor (1 — 2wi). A simultaneous fit to the mixing frequency and 
its amplitude allows the determination of both Amd and the mistag rates, Wi. 

Neglecting any background contributions, the probability density functions (PDFs) for the 
mixed (-|-) and unmixed (— ) events, TC±^sig, can be expressed as the convolution of the decay 



distribution for the i tagging category 

g-lAtl/r^O 

h±{At] Am^, Wi) = — [1 =F (1 — 2wi) cos A?n,rfAt] , 

Atbq 

with a time resolution function TZ{At — Attme', o,i), 

H±^sig{At] Am^, Wi, di) = h±{Attrue] Arud, Wi) ® TZ{At - Attrue] di). 

where Oj are the parameters of the resolution function. A log-likelihood function is then constructed 
from the sum of \og7i±^sig over all mixed and unmixed events, and over the different tagging 
categories. 

The log-likelihood is maximized to extract the mixing parameter Arrid and, simultaneously, 
the mistag rates, Wi. The correlation among these parameters is small, because the rate of mixed 
events at low values of At, where the B^B^ mixing probability is small, is principally governed by 
the mistag rate. Conversely, the sensitivity to Arud increases at larger values of At; for At near 
twice the lifetime, half of the B^ mesons will have oscillated. 

3 The BaBAR detector 

The BABAR detector is a 47r charged and neutral spectrometer described in more detail elsewhere Q. 
Charged particles are detected and their momenta measured by a combination of a 40-layer drift 
chamber (DCH) and a five-layer silicon vertex tracker (SVT) embedded in a 1.5-T solenoidal mag- 
netic field. Decay vertices are typically reconstructed with a resolution along the boost direction 
of 65 /im for fully reconstructed B mesons. A ring imaging Cherenkov detector, the DIRC, is used 
for charged hadron identification. A finely segmented CsI(Tl) electromagnetic calorimeter (EMC) 
is used to detect photons and neutral hadrons, and also for electron identification. The iron flux 
return (IFR) is segmented and instrumented with multiple layers of resistive plate chambers for 
the identification of muons and long-lived neutral hadrons. 

4 Event selection and B reconstruction 

The analysis uses a sample of multihadron events, which are selected by demanding a minimum of 
three reconstructed charged tracks and a total charged and neutral energy greater than 4.5 GeV in 
the fiducial region of the detector, and a reconstructed event vertex within 0.5 cm of the measured 
interaction point [^ in the plane transverse to the beamline. 

Electron candidates must satisfy a cut on the ratio of calorimeter energy to track momentum 
of 0.88 < E/p < 1.3, a cluster shape consistent with an electromagnetic shower, and DCH dE/dx 
and DIRC Cherenkov angle consistent with an electron. 

Muon candidates must satisfy requirements on the number of interaction lengths of IFR iron 
penetrated of Nx > 2.2, on the difference in the measured and expected interaction lengths pene- 
trated of N^^-'' — Nx < 1, on the position match between the extrapolated DCH track and the IFR 
hits, and on the average and spread of the number of IFR hits per layer. 

Pairs of photons in the EMC with energy above 30 MeV are constrained to the known vr*^ mass if 
they are within ±20 MeV/ c^ of the nominal invariant mass y], and their summed energy is greater 
than 200 MeV. 



K^ -^ Tr~^TT^ candidates are required to have an invariant mass between 462 and 534MeV/c^, 
and a x^ probability for the vertex fit of greater than 0.1%. The transverse flight distance of the 
Kg candidate from the primary event vertex must be greater than 2 mm. 

D^ candidates are identified in the decays channels K~^7r~ , K~^TT~7r^, i^'^vr^vr^vr" and i^^vr+vr". 
D^ candidates are selected in the K~^TT~7r~ and KgTT~ modes. Kaons from D^ decays and charged 
daughters from D^ — > K~^ tt~ are required to have a momentum greater than 200MeV/c. All 
other charged D daughters are required to have a momentum greater than 150 MeV/c. For D^ -^ 
K'^TT~iT^, we only reconstruct the dominant resonant mode D^ -^ K~^p^ , followed by p~ -^ it~-k^. 
The 7:~iT^ mass is required to lie within ibl50MeV/c^ of the nominal p mass Q and the angle 
between the vr" and D^ in the p rest frame, ^^o^i must satisfy |cos^^o^| > 0.4. D^ and D~ 
candidates are required to have momentum greater than 1.3 GeV/c in the T{4S) frame, an invariant 
mass within ±3a of the nominal value H and a x^ probability of the topological vertex fit greater 
than 0.1%. A mass constraint is applied to selected D^ candidates. 

D*~ candidates are formed by combining a D^ and a pion with momentum less than 450 MeV/c. 
The soft pion is constrained to originate from the beamspot when the D*~ vertex is computed. 
Those D*^ candidates with m{D^Tr~) — m{D^) lying within ±2.5(7 of the nominal value Q are 
selected, where a = l.lMeV/c^ for D^ — > K^it~tt^ mode and 0.8MeV/c^ for all other modes. 

J/V' — > e~^e~ or p~^p~ candidates must have at least one decay product positively identified 
as an electron or muon. Electron candidates outside the calorimeter acceptance must have DCH 
dE/dx information consistent with that for an electron. The second muon candidate, if within the 
acceptance of the calorimeter, must be consistent with being a minimum ionizing particle. J/ip 
candidates are required to lie in the invariant mass interval 2.95 (3.06) to 3.14GeV/c^ for the e~^e~ 
{p^p~) channel. 

B^ candidates in the flavor eigenstate decay modes D^*''7r^/p^/af are formed by combining a 
D*^ or D~ candidate with a 7r+, p+ (/9+ -^ 7r+7r°) or a^ {a'l -^ 7r+7r"7r+); likewise B^ -^ J/ipK*'^ 
candidates are reconstructed from combinations of J/f/' candidates with a K*^ (K*^ — > K~^7t~). 

For B^ -^ D*~p~^, the vr" from the p^ decay is required to have an energy greater than 300 MeV. 
For B^ -^ D*^al , the a^ is reconstructed by combining three charged pions, with invariant mass 
in the range 1.0 to 1.6GeV/c^ and a x^ probability of the vertex fit of the a'\^ candidate of greater 
than 0.1%. For most B^ modes, no particle identification or only a loose requirement is enough to 
achieve reasonable signal purities. 

Continuum background is rejected by requiring the normalized second Fox- Wolfram moment Q 
be less than 0.5. Further suppression is achieved by a mode-dependent restriction on the angle, 0th, 
between the thrust axes of decay products from B^ec and Btag respectively in the T(4S') frame. 

B^ candidates are identified with the difference /S.E between the energy of the candidate and 
the beam energy E^^^ in the center-of-mass frame and the beam-energy substituted mass rriEs = 
\/(^bram)^ - {pTY- Those candidates with mES > 5.2 MeV/ c^ and A£; = within ±2.5 standard 
deviations (typically |A£'| < 40 MeV) are selected. If there is more than one candidate satisfying 
these conditions only the one with the smallest /S.E is retained. Finally, a topological vertex fit of 
the candidate must converge. 

5 Decay time difference determination 

The decay time difference. At, between B decays is determined from the measured separation 
along the z axis between the reconstructed B = Byec and flavor-tagging decay B = Btag vertices 
Az = Zrec — ztag- This Az is then converted into At using the known T(4S') boost and correcting 
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Figure 1: Distribution of tties for all B^ hadronic candidates decaying into flavor eigenstate modes. 



on an event-by-event basis for the direction of the B mesons with respect to the z direction in 
the T{4:S) frame. The resolution of the At measurement is dominated by the z resolution of the 
tagging vertex. The Btag decay vertex uses all tracks in the event except those incorporated in 
B^ec- An additional constraint is provided by including a calculated B^ag production point and 
three-momentum, determined from the three momentum of the Brec candidate, its decay vertex, 
and the average position of the interaction point and the T{4S) boost. Reconstructed K'^ or A 
candidates are used as input to the fit in place of their daughters in order to reduce bias due to 
long-lived particles. Tracks with a large contribution to the x^ ^^^ iteratively removed from the fit, 
until all remaining tracks have a reasonable fit probability or all tracks are removed. Only those 
candidates with \Az\ < 3.0mm and uaz < 400 /xm are retained. The distribution of m^o, for the 
surviving candidates is shown in Fig. ffl along with a fit with a Gaussian distribution for the signal 
and the ARGUS function [^ for the background. At this point the sample contains 6662 it 93 
signal events, with a purity, defined in the region ttzes > 5.27GeV/c^, that varies between 80-95% 
depending on the B^ decay mode. 

In the likelihood, the time resolution function can be approximated by a sum of three Gaussian 
distributions with different means and widths, 

3 , 
fc=i ^fc^^vr 
where, for the core and tail Gaussians, the widths ai^2 = 5'i^2 x cta* are the scaled event-by-event 
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measurement error, cr/\t: derived from the vertex fits. The third Gaussian, with a fixed width of 
C3 = 8ps, accounts for less than 1% of outher events with incorrectly reconstructed vertices. A 
separate core bias, 6i^i, is allowed for each tagging category i to account for small shifts due to 
inclusion of charm decay products in the tag vertex, while a common bias 82 is used for the tail 
component. The tail and outlier fractions and the scale factors are assumed to be the same for all 
decay modes, since the precision of the B^ag vertex dominates At. This assumption is confirmed 
by Monte Carlo simulation studies. 

6 Flavor tagging 

After the daughter tracks of the reconstructed B are removed, the remaining tracks are analyzed 
to determine the flavor of the Btag, and this ensemble is assigned a tag flavor, either B^ or B^. For 
this purpose, we use the flavor tagging information carried by primary leptons from semileptonic B 
decays, charged kaons, soft pions from D* decays, and more generally by high momentum charged 
particles to uniquely assign each event to a tagging category. The effective tagging efficiency 
Qi = £4(1 — 2wi)'^, where £i is the fraction of events assigned to category i, is used as the basis for 
optimization of category selection criteria. The statistical error on Am^ is proportional to Xj^fQ, 
where Q = Y^Qi. 

Events are assigned to the Lepton category if they contain an identified lepton with a center- 
of-mass momentum greater than 1.0 or l.lGeV/c for electrons and muons, respectively. The mo- 
mentum requirement selects mostly primary leptons by suppressing opposite-sign leptons from 
semileptonic charm decays. 

Kaons are identified with a neural network based on the likelihood ratios in the SVT and 
DCH, derived from dE/dx measurements, and in the DIRC, calculated by comparing individual 
photomultiplier hits with the expected pattern of Cherenkov light for either kaons or pions. The 
charges of all identified kaons are summed, and if ^ Qk 7^ the event is assigned to the Kaon 
category. 

The final two categories involve a multivariable analysis using a neural network, which is trained 
to identify primary leptons, kaons, and soft pions, and the momentum and charge of the track with 
the maximum center-of-mass momentum. Depending on the value of the output variable from the 
neural net, events are given a B^ or B^ tag and assigned to the mutually exclusive categories NTl 
(more certain tags) or NT2 (less certain tags). About 30% of all events are assigned to no tagging 
category and are excluded from the analysis. 

Tagging assignments for events are made mutually exclusive by the hierachical use of the tagging 
categories. Events with a Lepton tag and no conflicting Kaon tag use the Lepton category. If no 
Lepton tag exists, then the Kaon category is used, if a tag exists. Otherwise the two neural network 
categories are used. The number of tagged events per category is given in Table |l[ In total, there 
are 4538 it 75 tagged signal events. 

7 Background PDF 

In the presence of backgrounds, the probability distribution functions 'H±^sig must be extended to 
include a term for each significant background source, which are allowed to differ for each tagging 
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Table 1: Event yields for the different tagging categories used in this analysis, as obtained from 
fits to the m-Es distributions after all selection requirements. The purity is quoted for rriES > 
5.270 MeV/c^. 

Category Tagged Efficiency (%) Purity (%) 



Lepton 


754 ± 28 


11.3 ±0.4 


97.1 ±0.6 


Kaon 


2317 ± 54 


34.8 ± 0.6 


85.2 ±0.8 


NTl 


556 ± 26 


8.3 ±0.3 


88.7 ±1.5 


NT2 


910 ±36 


13.7 ±0.4 


83.0 ± 1.3 



Total 4538 ±75 68.1 ±0.9 86.7 ±0.5 



category: 

T~L±,i = fi,sigH±,sig + 2_^ fi,kB±,i,k{^t]b±^i^k) 
k=bkgd 

where the background PDFs, B±^i^k, provide an empirical description for the At distribution of the 
background events in the sample. The fraction of background events for each source and tagging 
category is given by /j^^, while 6±,i,fc are parameters used to characterize each source of background 
by tagging category for mixed and unmixed events. 

The probability that a B^ candidate is a signal or a background event is determined from 
a separate fit to the observed tit-es distributions of B^^c candidates in each of the four tagging 
categories. We describe the m-Es shape with a single Gaussian distribution 5(?tt,es) for the signal 
and an ARGUS parameterization ^(tties) for the background. Based on this fit, the event-by-event 
signal probabilities fsig^i are given by 



f ( N 5(mEs) 



SiruEs) ±^(mEs) 



The sum of signal and background fractions is forced to unity. 

The At distributions of the combinatorial background are assumed to be described with a zero 
lifetime component and a non-oscillatory component with non-zero lifetime. We fit for separate 
resolution function parameters for the signal and the background in order to minimize correlations 
of the time structure between background and signal. Candidates with low signal probability, i.e., 
in the ttt-es sideband region below 5.27GeV/c^, dominate the determination of these background 
parameters. 

8 Extraction of Arud 

The value of Arud is extracted from the tagged flavor-eigenstate B^ sample with an unbinned 
maximum likelihood fit involving a total of 34 parameters, including Am^. The value of Arud was 
kept hidden throughout the analysis until the systematic errors were finalized, in order to eliminate 
possible experimenter's bias. 

The results from the likelihood fit to the tagged B^ sample are summarized in Table 0. The 
probability of obtaining a likelihood smaller than the one observed is determined to be 28% from a 
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Table 2: Results from the likelihood fit to the At distributions of the tagged hadronic B^ decays. 
Arud and the mistag rates include small corrections corresponding to the difference between the 
generated and reconstructed values in simulated signal events. 



Parameter 
Arrid (hps^^) 

''^Lepton 
''WKaon 
■W^NTl 
U'MT2 

Si 

S2 



Value 



Parameter 



Value 



0.519 ±0.020 
0.085 ±0.018 
0.167 ±0.014 
0.195 ±0.026 
0.326 ± 0.024 
1.42^ 
5.5^ 



hO.08 
-0.09 
h2.0 
-1.6 



'^l,Lepton(ps) 

'^l,Kaon(ps) 

'^l,NTl(ps) 

'^1,NT2(PS) 

52(PS) 

/2 

/3 



+0.11 ±0.07 
-0.20 ±0.05 
±0.01 ± 0.09 
-0.20 ±0.08 
-1.2^ 
0.032^ 
0.001 ± 0.004 



hO.9 
-1.6 
-0.03 
-0.02 
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Figure 2: At distributions in data for the selected mixed and unmixed tagged hadronic B^ decays 
("^ES > 5.27GeV/c^), with overlaid the projection of the likelihood fit. The background contribu- 
tion is indicated by the dashed curve. 
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Figure 3: The time-dependent asymmetry Amixi^i) between unmixed and mixed events as a 
function of lAtl. 



large number of simulated experiments, each generated according to the parameters obtained from 
the fit. The measured value for Amd is: 

Arud = 0.519 ± 0.020(stot) ± 0.016(syst) /ips"^ 

where the sources of systematic error are discussed below. The observed distribution of mixed and 
unmixed events and the asymmetry, Amix are shown in Fig. |^ and ^ as a function of At along with 
projections of the likelihood fit result. 

9 Systematic studies and cross checks 

The conversion of Az to At introduces a systematic uncertainty (±0.007 /ips~^) due to the limited 
knowledge of the PEP-II boost, the z length scale of BABAR (determined from secondary interactions 
in a beam pipe section of known length) and the Bj-^c momentum vector in the T(45) frame. 

The signal At resolution parameters are determined directly from data by the fit, contributing 
±0.006 in quadrature to the statistical error. Residual uncertainties (±0.005 h ps~^), are attributed 
to the choice of the parameterization, its description of the outliers, and the capability of the 
resolution model to deal with various plausible misalignment scenarios applied to the Monte Carlo 
simulation. 

The parameters of the background At distribution are left free in the likelihood fit, but sys- 
tematic errors (±0.005 /ips~^), are introduced by the residual uncertainty from the r^ES fit used 
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to determine the signal probability, the assumed parameterization of the background At distribu- 
tions and resolution function, and the small amount of correlated B^ background remaining in the 
sample. 

Finally, statistical limitations of Monte Carlo validation tests (±0.004 ^ps~^), the full size of 
a (negative) correction obtained from Monte Carlo (±0.009 ?ips~^), and the variation of the fixed 
B^ lifetime (due to its negative correlation with Amd) within known errors Q (±0.006 7ips~^) 
contribute. A summary of these sources of systematic error for the hadronic B^ sample is shown 
in Tables |. 

Various checks on the consistency of the result were performed by splitting the data in sub- 
samples according to several key variables, including (but not limited to) the decay modes of Sreo 
the tagging category of Btag , and the flavor of either Brec or Bfag ■ The value of Amd was found 
to be consistent for all sub-samples. 

Table 3: Systematic uncertainties for Amd- 



Source 


SAnid \h ps -^] 


Beamspot 


position and size 


0.001 


Az to At conversion 




PEP-II boost 


0.003 


z scale 


< 0.005 


method 


0.004 


At resolution 




outliers 


0.002 


parameterization 


0.003 


SVT alignment 


0.004 


Backgrounds 




At model 


0.001 


Resolution parameterization 


0.003 


Fractions 


0.003 


Correlated 


0.002 


Monte Carlo 




statistics 


0.004 


correction 


0.009 


B^ lifetime 


0.006 


Total Systematic Error 


0.016 



10 Summary 

We have measured the value for Arud to be: 

Amd = 0.519 ± 0.020(stot) ± 0.016(syst) /ips'^ 

This result is one of the best single measurements available and is consistent with the current world 
average Q. Moreover, the error on Anid is still dominated by the statistical size of the reconstructed 
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B^ sample, leaving substantial room for further improvement as more data is accumulated at 
BABAR. The measurement shares the same flavor-eigenstate sample as used for the determination 
of sin2/3. Thus, it provides an essential validation for the reported sin2/3 result H and, in particular, 
the mistag rates that appear as coefficients of the mixing asymmetry. 
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